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The antiferromagnetic Ising model on a checkerboard lattice has an ice-like ground state manifold 
with extensive degeneracy, and, to leading order in J xy , deconfined spinon excitations. We explore 
the role of cyclic exchange arising at order J^y/Jz on the ice states and their associated spinon 
excitations. By mapping the original problem onto an equivalent quantum six-vertex model, we 
identify three different phases as a function of the chemical potential for flippable plaquettes — a 
phase with long range Neel order and confined spinon excitations, a non-magnetic state of resonating 
square plaquettes, and a quasi-collinear phase with gapped but deconfined spinon excitations. The 
relevance of the results to the square-lattice quantum dimer model is also discussed. 

PACS numbers: 75.10.-b, 75.10.Jm, 64.60.Cn 



The past decade has seen a great renaissance in the 
study of frustrated quantum spin systems. On the ex- 
perimental front, advances in the synthesis of magnetic 
oxides have given rise to a great wealth of new frustrated 
materials with highly unusual and interesting properties. 
And, at the same time, highly frustrated models have be- 
come a favorite playground of theorists seeking to under- 
stand unconventional phase transitions and excitations. 

Recently, it was proposed that the geometric frustra- 
tion present on the pyrochlore lattice could give rise to 
fractional charges in two or three dimensions^ in a phys- 
ically realistic model based on strong nearest neighbor 
repulsion close to commensurate filling >S The charge or- 
dering problem considered in .1] is classically equivalent 
to one of Ising antiferromagnetism, and in this paper we 
consider the simplest possible test case for these ideas, 
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FIG. 1: a) The checkerboard lattice on which the ice states, 
and b) the square lattice on which the states of the six-vertex 
model are defined. Any Ising state obeying the ice rules, 
e.g. that shown in (a), is equivalent to (b) six-vertex model 
configuration. In the state shown, the upper left corner has 
Neel order, while the lower right corner has collinear order. 
Flippable plaquettes are denoted with circles. In the case of 
the six-vertex model, these have a definite sense of rotation 



the XXZ Heisenberg model on a checkerboard (2D py- 
rochlore) lattice. We proceed by mapping this model 
onto an equivalent, quantum six-vertex model (Q6VM), 
and describe the nature of the ground state and low lying 
spin excitations of this model as a function of a control 
parameter V, which acts as a chemical potential for those 
"flippable plaquettes" accessible to cyclic exchange. 

We identify three different ground states, a phase with 
long range Neel order, a non-magnetic state of resonat- 
ing square plaquettes, and a partially disordered phase of 
"isolated states" with extremely large ground state de- 
generacy, referred to as the "quasi-collinear" phase be- 
low. Because of the anisotropy of the model, all spin 
excitations are gapped. It is possible to identify the low- 
est lying excitations of the Neel phase as spin waves, and 
those of the quasi-collinear phase as deconfined spinons. 
We also identify the special role of the isolated states 
in supporting fractional excitations. Many of these re- 
sults are also relevant to the much studied square lattice 
quantum dimer model (QDM)i^ 

Model and mapping onto Q6VM. We take as a start- 
ing point the spin-1/2 anisotropic Heisenberg model with 
antiferromagnetic interactions, J z ,J xy > 0, in the limit 

Jz > J'. 
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Here the sum run s over the bonds of the 2D py- 

rochlore or checkerboard lattice, shown in Fig.^Jt). In the 
Ising limit, J xy — 0, this model has an extensive ground 
state degeneracy — every state with exactly two up and 
two down spins per tetrahedron (cross linked square) is a 
ground state. For historical reasons, this is known as the 
"ice rules" constraint. Topologically, "ice" states have 
the structure of closely packed loops of up and down 
spins, and are separated by a gap J z from the lowest 
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lying excited state. Flipping any given down spin con- 
nects two adjacent loops of up spins, creating two "T- 
junction" like topological defects (spinons), which prop- 
agate independentlyi^ The pyrochlore (checkerboard) 
lattice is bipartite in tctrahedra. Spinons are created 
in A and B sublattice pairs, and move so as to preserve 
tetrahedron sublattice. 

By drawing an arrow from the center of A to B sublat- 
tice tetrahedra where they share an up spin, and from B 
to A where they share a down spin, one can show that the 
many ground states of the Ising model on a checkerboard 
lattice are in exact, one-to-one correspondence with the 
states of the classical six vertex model (6VM)j^*£ widely 
studied as a 2D analogue of water ice. From this map- 
ping, we know that a) the ground state manifold of the 
Ising model grows as W cc (4/3) 3Ar / 4 where N is the 
number lattice sitesi and b) all correlation functions de- 
cay algebraically, 8 

Up to this point, our analysis contains only classical 
statistical mechanics and simple topological arguments. 
Quantum mechanics reenters the problem when we con- 
sider a small but finite J xy C J 2 - In this case, the ice 
states are no longer eigenstates. Short lived virtual exci- 
tations enable the system to tunnel between different ice 
state configurations wherever pairs of up spins and down 
spins occur diagonally opposite one another on one of 
the empty square plaquettes of the checkerboard latticed 
The allowed reconfigurations of these "flippable plaque- 
ttes" can be described within degenerate perturbation 
theory by the effective Hamiltonian 

Ti-2nd = — y~ y](StS 2 S^S 4 + S t S^Sz S£) , (2) 
Jz □ 

where the indices 1 to 4 count consecutive sites (either 
clockwise or anticlockwise), of an empty plaquette^ 

In terms of the 6VM representation, Eq. J5J acts on a 
plaquette where four arrow are joined nose to tail, so as to 
invert all of the arrows and change the sense of rotation 
of the plaquette (c.f. [To|'). The quantum dynamics in 
the Q6VM we consider are directly analogous to the res- 
onance of dimers in the QDM^ studied as an simplified 
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FIG. 2: The phase diagram of the model as a function of 
V/t. The Neel phase breaks the point group, while the plaque- 
tte phase breaks the translational symmetry. The Rokhsar- 
Kivelson point is marked RK. 
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FIG. 3: Energy level diagram of the 32-site pyrochlore-slab 
with periodic boundary conditions as a function of V/t, ob- 
tained by numerical diagonalization. We have shown the 
first 8 levels. Inset: the first two excited states cross at 
V c /t = —0.3437 (the axes are the same as of the main plot). 

model of a resonating valence bond stated Formally, in 
fact, the Hamiltonian is exactly the same, although the 
Hilbert space on which it acts is different. And, as in the 
QDM, we anticipate that quantum effects will in general 
select a ground state with finite degeneracy from the vast 
manifold of classically-allowed ice states. 

As such, there is only one (kinetic) energy scale in the 
problem, t = J xy /J z . However in order to study the 
different possible phases of the model it is useful to in- 
troduce a further control parameter. A suitable control 
parameter for the QDM is a diagonal term which counts 
the number of dimers which can resonate in any given 
dimer covering. By direct analogy, we introduce a diago- 
nal interaction V to the Q6VM which counts the number 
of flippable plaquettes 

n = J2[v(\o)(o\ + \cj){o\)-t(\o){cj\ + \cj)(o\)} , 
□ 

(3) 

where the \0) and |0) states represent squares with 
the respective circular arrow configuration on the square 
edges, as seen in Fig. ^b). We note that, for a system 
with periodic boundary conditions, the net flux of ver- 
tex arrows through any given horizontal or vertical cut 
defines a set of winding numbers which are conserved by 
the Hamiltonian © . 

Our approach to determining the different phases of 
the Hamiltonian J2J is the numerical diagonalization of 
clusters with periodic boundary conditions of up to 64 
spins, within the ice rules manifold of states, supple- 
mented with topological and symmetry arguments. De- 
tails of these, together with further analysis of the related 
fermionic charge-ordering problem will be discussed fur- 
ther in separate publicationsiiSii^ 

Phase diagram. We first consider the nature of the 
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ground state as a function the chemical potential for flip- 
pable plaquettes, V. Our results are summarized in the 
phase diagram Fig. and the numerical evidence for 
each phase discussed below. 

Negative values of V favor states with flippable plaque- 
ttes. The state with the greatest possible number of flip- 
pable plaquettes is the Neel state, and this must be the 
grounds state for V — > — oo. The Neel state is two fold 
degenerate in the thermodynamic limit. For finite V/t, in 
a finite size system, quantum fluctuations lift this degen- 
eracy, as seen in the low-energy spectrum of the Q6VM 
(Fig. El). We find a single phase for V < — 0.3i, which 
we identify as the Neel phase. Both the the symmetric 
and antisymmetric combinations of the two symmetry- 
breaking Neel ground states are visible in the spectrum, 
marked "GS" and "Neel" respectively. At a value of 
V ~ — 0.3t, a third energy level, marked "Plaq" crosses 
the first excitation "Neel" . We interpret this as evidence 
for a quantum phase transition into a resonating plaque- 
tte phase, discussed below. From finite size scaling of 
the spectrum (Fig.^J we estimate the critical value to be 
V c = — 0.3727i in the thermodynamic limit. As the com- 
peting Neel and plaquette order parameters break lattice 
symmetries in very different ways, the transition between 
them is presumably of first order. 

We find a single phase extending from — 0.3i < V < t, 
including the XXZ point V = 0. This phase terminates 
in the special high symmetry point V = t for which the 
Hamiltonian © of the Q6VM can be written as a sum 
of projection operators : 

n RK = tJ2(\o)-\o))((o\-(o\) . (4) 
□ 

Following Rokhsar and Kivelson (RK) ^ we can construct 
a zero eigenvalue state of the 7Yrk by taking the linear 
combination of all the states in a given topological sector 
with the same amplitude. Since this state is annihilated 
by the positive semi-definite Hrk, it must be a ground 
state. As in the QDM, static correlations can be com- 
puted exactly at this point. Like the correlation functions 
of the 6VM, they decay algebraically with distance. 

At the RK point, kinetic and potential energy are per- 
fectly balanced; in the plaquette phase kinetic energy 
dominates, and resonating plaquettes repel one another 
so as gain the maximum kinetic energy^ The result- 
ing state is essentially a Peierls-like distortion of the RK 
state in which only A (B) sublattice plaquettes resonate 
nn 1(B) (|0)+ |0)). The way in which the phase breaks 
lattice symmetries — it is two-fold degenerate, and in- 
variant under operations which map the alternating A 
and B plaquette sublattice onto themselves — suggest 
the plaquette phase of the Q6VM is an Ising analogue 
of the a SU(2) valence-bond crystal of resonating pla- 
quettes. Such a phase has been proposed in the context 
of the square lattice QDM^ Furthermore, the ground 
state of the Heisenberg-model on a checkerboard lattice 
is a valence bond crystal of SU(2) singlets formed on al- 
ternate empty plaquettes ji£ with a possibility of an adi- 
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FIG. 4: Estimate of the phase boundary between the Neel 
and plaquette phases. Empirically, the values of V c where the 
level crossings occur scale as —0.3727 + 1.86-/V~ 3//2 . Values 
are shown for 32, 36, 40, 52, and 64 pyrochlore-slab sites. 



abatic continuity between the ground state of the XXZ 
and SU(2) symmetric Heisenberg models^ 

For V > t the ground state is the highly degenerate 
manifold of "isolated" states with no flippable plaquettes. 
They are eigenstates with energy for any value of V/t, 
and become the ground state for V > ti& The prototype 
of an isolated state is the collinear configuration shown 
in Fig ^ In this reference state all vertex arrows point 
from left to right or from top to bottom. Inverting the di- 
rection of the arrows along an arbitrary number of lines, 
subject to the constraint that all of them are either hori- 
zontal or vertical, creates new isolated states. This leads 
to a ground state degeneracy which grows as 4(2 P — 1) 
for regularly shaped clusters, where p ~ y/N. In these 
states, the direction of arrows along either the horizon- 
tal or vertical lines is long-range ordered, but quantum 
effects none the less fail to select a ground state with fi- 
nite degeneracy. We refer to this phase of the Q6VM as 
"quasi-collinear" . Finally, since the transition between 
the quasi-collinear phase and the resonating plaquette 
phase takes place through the softening of specific exci- 
tation (discussed below), we identify it as second order. 

Excitations. First let us first consider the nature of 
excitations at fixed S z — 0. A state with n flippable 
plaquettes has a diagonal matrix element nV and is con- 
nected to n other states. Gerschgorin's theorem places 
a bound \Hu — e»| < J2j l-^ijl on the separation of the 
i-th eigenvalue Si from the diagonal matrix element Ha . 
In the case in point, this bound is \nV — e»| < nt, or 
n(V — t) < £i < n(V + 1). The smallest energy in an ar- 
bitrary topological sector is thus larger than V — t, which 
gives a lower bound on the value of the gap in the quasi- 
collinear phase for V > t. This above argument permits 
a gapless spectrum at the RK point V = t. In fact it is 
possible to explicitly construct a family of states with a 
gap that vanishes at the RK point, as shown in Fig. 0a): 
the energy spectrum of this particular excitation forms a 
continuum between 2V — 2t and 2V + 2t. 

Now let us consider spin excitations with S z — ±1. 
If we neglect virtual processes at order J% y /J z , and the 
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FIG. 5: (a) A "leap-frog" excitation in the quasi-collinear 
phase. Two flippable plaquettes (denoted with circles) are 
created by reversing the arrows of the collinear reference state 
on a line with a single one-step kink. The motion of the pair of 
flippable plaquettes is equivalent to a one-dimensional hop- 
ping model with an energy spectrum e(k) = 2V + 2tcosfc, 
where k is an effective one-dimensional momentum, (b) and 
(c): The deconfined spinons in the collinear phase (black 
dots). Note that spinons hop so as to stay within a given 
sublattice of the (bi-partite) square lattice. 



possibility of entropic confinement at finite temperature, 
these propagate as independent fractional excitations. 1 
Quantum effects may, or may not, act to confine these 
excitations, depending on the type of correlations present 
in the ground state they select. The Neel ground state 
has a two fold ground state degeneracy, and separating 
the topological defects created by flipping a spin creates 
a string of unflippable plaquettes. This leads to confine- 
ment of spinons, and the low lying spin excitations of the 
Neel phase of our model have the same quantum num- 
bers as a spin wave. On general grounds, we expect the 
same to be true of the plaquette phase. 

The manifold of isolated states selected by V can sup- 



port deconfined spinons, however. Since no new flippable 
plaquettes are introduced into isolated states by flipping 
a single spin, and the pair of topological defects created 
by flipping a single spin can be separated without cre- 
ating new flippable plaquettes, spinons are deconfined. 
An example of a pair of deconfined spinon excitations is 
shown in Fig.[5Jb)-(c). For V 3> J xy ^> t, spinon motion 
is movement is confined to the x and y directions, but by 
scattering off one another, a pair of spinons can explore 
the full two dimensional space of the lattice. Whether 
a more general class of deconfined spinon excitation be- 
comes possible as one approaches the RK point remains 
an open question. We also note that while fixing the 
boundary conditions will lift the degeneracy of the iso- 
lated state manifold, it need not affect the arguments for 
spinon deconfinement presented above. 

Conslusions We have established that, as a function of 
the chemical potential for "flippable plaquettes" accessi- 
ble to cyclic exchange, the XXZ Heisenberg model on a 
checkerboard lattice exhibits Neel, resonating plaquette 
and quasi-collinear phases. If virtual processes at J% y / J z 
are ignored, spinon excitations in the XXZ Heisenberg 
model are deconfined. We have shown explicitly that a 
subset of spinon excitations — those associated with iso- 
lated states — remain deconfined even when these quan- 
tum effects are taken into account. Finally, we mention 
that the equivalents of both the "leap-frog" and spinon 
excitations can also be constructed in the square-lattice 
QDM for V > 
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